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We  have  investigated  transient  heat  transfer  to  liquid  hydrogen.  Thin  carbon  films 
and  Pt  foils  submerged  in  liquid  hydrogen  received  stepped  inputs  of  power  ranging 
from  1  to  42  V/cm2,  and  the  onset  of  nucleate  or  film  boiling  was  obtained  for  each 
power  level.  The  critical  heat  flux  was  found  to  be  approximately  8  V/cm2,  with  the 
transition  to  film  boiling  occurring  in  times  less  than  10*3  s.  Premature  film  boiling 
can  be  related  to  the  positive  temperature  coefficient  of  resistance  and  the 
narrowness  of  the  heaters.  Thermometric  devices  and  power  generation  equipment 
selection  are  discussed. 


1.  INTRODUCTION 

The  goal  of  this  study  is  to  fill  the  void  in  data  for  transient  heat  transfer  in 
liquid  hydrogen.  The  experimental  technique  has  been  previously  used  in  liquid  helium 
[1,2,3]  and  liquid  nitrogen,  [1,4,5]  and  uses  stepped  inputs  of  power  to  submerged 
heaters.  The  duration  of  the  input  ranges  from  10*5  to  10  s,  with  the  onset  of  nucleate 
boiling  occurring  in  times  on  the  order  of  10'3  s.  The  surface  temperature  is  monitored 
throughout  the  period  of  application  of  power. 

In  Section  2  the  report  describes  the  selection  and  assembly  of  the  experimental 
apparatus,  including  the  safety  considerations  for  using  liquid  hydrogen  at  NIST,  and 
the  various  heater/theraometers  used.  An  analysis  of  transient  heat  transfer  behavior 
is  presented  in  Section  3,  followed  by  a  presentation  of  experimental  results  in 
Section  4.  Finally,  the  conclusions  are  given  in  Section  5. 


Keywords:  boiling;  conduction;  hydrogen;  premature  film  boiling;  thermometry; 

transient  heat  transfer 
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2.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Procedures  used  in  these  experiments  were  subject  to  the  NIST  guidelines  for  the 
laboratory  use  of  liquid  hydrogen  [6] .  The  prevention  of  backflow  of  air  into  the 
apparatus,  devar  breakage,  and  the  normal  venting  from  the  experiment  were 
specifically  identified  as  safety  concerns.  Because  this  study  focused  on  transient 
behavior,  many  potential  hazards  were  minimized  by  the  size  of  the  experiment.  Less 
than  2  liters  of  liquid  hydrogen  were  used  in  any  one  experiment ,  and  the  duration  of  an 
experiment  was  less  than  hair  a  day.  Also,  experiments  in  liquid  nitrogen  enabled  us 
to  identify  potentially  hazardous  techniques  which  were  then  modified  for  use  with 
liquid  hydrogen. 

The  experimental  technique  used  for  this  study  was  the  application  of  a  step  in 
power  applied  to  a  thin  film  resistance  device.  The  change  in  resistance  of  this 
device  was  used  to  determine  the  temperature  change  and  heat  transfer  characteristics 
at  the  its  surface  and  .  The  selection  of  the  appropriate  electronic  equipment  and 
test  devices  is  summarized  in  this  section. 

2.1  Overview  of  the  Experimental  Apparatus 

A  small,  vacuum—  jacketed  dewar  with  optical  viewports  was  used  as  the  cryostat. 
Liquid  nitrogen  shielded  the  dewar.  The  test  section  was  mounted  vertically  at  the  end 
of  a  long  rod  in  the  neck  of  the  dewar.  The  temperature  of  the  bath  was  measured  with  a 
germanium  resistance  thermometer.  The  assembly  is  shown  in  Figure  1.  The  dewar  was 
surrounded  with  a  plexiglass  shield  and  metal  curb,  which  served  as  protection  from 
glass  fragments  from  dewar  breakage.  0-ring  seals  were  used  on  the  flange  fitting  for 
the  main  dewar  compartment  and  for  the  electrical  leads . 

During  an  experiment  the  vent  gas  left  the  dewar  and  was  warmed  to  ambient 
temperature  within  a  large  copper  coil,  passed  through  two  check  valves,  and  exited 
above  the  rooftop.  The  check  valves  prevented  backflow  of  air  into  the  apparatus.  An 
electrochemical  hydrogen  detector  was  located  above  the  dewar.  The  instrument  rack 
was  purged  with  fresh  air  during  experiments. 

2.2  Power  Generation  Devices 

An  ideal  step  or  pulse  input  of  power  has  a  very  short  rise  time  with  a  minimal 
overshoot,  a  period  of  steady  power,  and  a  fast  shutoff  with  a  minimal  backswing.  Most 
power  supplies  are  designed  for  a  specific  resistive  load  in  order  to  produce  the  ideal 
waveform.  Changing  the  size  of  the  resistive  load  can  significantly  distort  the  shape 
of  the  power  input.  Consider,  too,  that  limitations  in  the  range  of  a  power  supply  can 
dictate  what  the  resistance  of  the  load  must  be.  Finding  acceptable  solutions  to  this 
trial-and— error  problem  is  most  important  in  obtaining  useful  data.  Several  power 
generators  vere  used  in  this  study  to  attempt  to  match  impedances  with  the  various  test 
sections.  Following  is  a  summary  of  the  various  power  supplies  used. 

Square  wave  pulse  generator.  This  pulse  generator  was  built  by  the  Fracture  and 
Deformation  Division  of  NIST.  It  had  a  pulse  width  of  either  1(H  or  10'5  s,  a  maximum 
current  of  20  A  and  voltage  output  of  50  V,  and  was  designed  to  handle  a  resistive  load 
of  2.5  to  2.6  ohms.  This  unit  was  limited  to  the  production  of  repetitive  pulses  which 
were  unfortunately  characterized  by  large  negative  backswings;  hence,  it  was  used  only 
to  obtain  preliminary  data. 
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High  power  Dulse  generator.  This  commercially  available  generator  could  produced 
either  single  or  repetitive  pulses  (up  to  30  kHz)  and  variable  pulse  widths 
(approximately  10*  to  10*  s),  by  using  different  plug-in  modules  to  set  the  output. 
The  voltage  ranged  from  48  to  1000  V  and  the  current, from  5  to  100  A,  depending  on  the 
module.  Tne  pulses  produced  by  this  apparatus  exhibited  an  overshoot  of  up  to  157,  and 
negative  backsving,  but  both  characteristics  were  modulated  by  proper  choice  of  load 
impedance  and  by  the  addition  of  several  electrical  components  on  the  output.  The 
electrical  circuit  is  shown  in  Figure  2.  The  largest  drawback  of  this  unit  was  the  duration 
of  the  pulse.  It  was  found  that  constant  power  could  not  be  maintained  for  longer  than 
10*  s,  which  is  the  minimum  time  required  to  reach  steady-state  boiling  in  either 
liquid  nitrogen  or  liquid  hydrogen.  Also,  some  of  the  output  units  exhibited  a  ramping 
output  power  atypical  of  specifications.  This  unit  was  limited  for  use  only  in 
producing  short,  high  power  pulses. 

Square— wave  pulse  generator  with  built-in  delay.  Ve  also  used  on  a  limited  basis  a 
third  pulse  generator  fabricated  in  our  Division  for  other  transient  heat  transfer 
experiments.  This  unit  has  a  single  1  to  100  A  square  pulse  of  10*s  to  1  s  duration. 
This  generator  is  incapable  of  producing  pulses  shorter  than  10*  s  or  less  than  1  A  in 
amplitude.  Since  it  was  an  integral  component  of  another  experiment,  its  availability 
was  restricted. 

Bipolar  power  amplifier  with  function  generator.  A  bipolar  power  amplifier  with  a 
function  generator  was  used  to  provide  step  inputs  of  low  to  medium  power.  Its  maximum 
voltage  is  72  V.  Power  inputs  of  0.01  to  55  V/cmJ  were  used.  This  unit  has  a  maximum 
current  output  of  5  A,  which  was  found  to  be  insufficient  for  some  tests.  Since  it  can 
operate  in  a  fast  mode,  with  ramp  times  of  less  than  2.5  x  10*  s,  it  was  able  to  produce 
clean,  square,  stepped  input  of  power. 

2.3  Transient  Response  Thermometry 

The  selection  of  thin  film  resistance  devices  for  transient  response  thermometry 
was  guided  by  several  factors.  It  is  desirable  to  select  a  material  which  has  a  fast 
response  as  a  heater  and  whose  resistance  is  a  known  function  of  temperature.  A  single 
device  used  as  both  the  heater  and  thermometer  would  eliminate  the  temperature 
difference  and  time  lags  in  heat  transfer  between  the  heater  and  thermometer. 
Platinum  was  identified  as  a  potential  material  because  of  its  suitable  resistance 
response  in  the  20  K  range  and  its  ease  of  fabrication  into  usable  devices.  Carbon 
films  demonstrated  characteristics  which  proved  beneficial  to  obtaining  fast 
temperature  responses.  For  all  devices  used,  significant  study  was  required  to  match 
the  impedance  of  the  heater/thermometer  with  that  of  the  power  generator,  while 
optimizing  the  heat  transfer  geometry.  This  section  will  discuss  our  experience  with 
each  type  of  heater/thermometer,  briefly  indicating  the  suitability  of  each  device. 

Vapor-deposited  films.  Several  vapor-deposited,  platinum  thin  film  resistance 
devices  on  quartz  substrates  were  fabricated  to  our  design  by  the  Cryoelectronics 
Division  at  NIST  in  Boulder  to  serve  as  heaters  and  as  thermometers.  A  schematic  of  a 
film  is  shown  in  Figure  3.  The  thickness  of  the  devices  ranged  from  1  to  4  x  10*  cm. 
The  resistance  was  adjusted  by  increasing  path  length  of  the  platinum  with  a 
serpentine  pattern.  These  test  sections  exhibited  poor  sensitivity  at  low 
temperatures,  with  an  average  resistivity  ratio  Rrt/RRm  at  77  K  of  0.6.  Although 
annealing  the  thin  films  at  550  °C  appeared  to  relieve  fabrication  stress  points,  the 
resistivity  ratio  at  77  K  improved  only  to  0.45.  The  resistivity  ratios  at  20  K  and  4  K 
were  both  0.25,  indicating  no  sensitivity  at  these  lower  temperatures.  This  behavior 
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Figure  2.  Schematic  of  electrical  circuit  for  transient  hydrogen  heat  transfer 
experiments . 


Figure  3.  Schematic  of  vapor— deposited  Pt  film  on  a  quartz  substrate  for  use  as  a 
neater  and  thermometer . 
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is  similar  to  that  of  an  impure  metal  and  suggests  that  the  mean  free  path  of  the 
electrons  is  larger  than  the  thickness  of  the  film.  Attempts  to  increase  the  thickness 
of  the  films  failed.  Consequently,  these  heater/thermometers  could  not  be  used  in 
these  experiments. 

Platinum  foil  on  ouartz  and  fiberglass-epoxv  substrates.  Pieces  of  platinum  foil  4  x 
10'4  cm  thick  were  mounted  on  quartz  or  fiberglass  epoxy  substrates  and  cut  into 
serpentine  patterns  to  reach  a  specific  resistance.  A  schematic  drawing  of  this 
heater/thermometer  is  shown  in  Figure  4.  The  resistivity  ratios  at  77  K  and  4  K  for 
unsupported  samples  of  foil  were  BRrt/RKm  =  0.185  and  0.0138,  respectively. 
Preliminary  data  for  heat  fluxes  less  than  20  V/cm2  in  saturated  liquid  hydrogen  showed 
fast  response  and  high  sensitivity. 

There  were  problems  with  this  design.  These  thermometers  had  to  be  "trained" 
before  they  would  hold  a  steady  resistance  at  low  temperatures.  After  cycling  20  to  30 
times  between  LNj  and  room  temperature,  the  resistance  would  stabilize  at  77  K  and 
below.  Ve  postulate  that  the  difference  in  the  thermal  expansion  between  the 
substrate  and  the  platinum  foil  induced  stress  points  in  the  serpentine  structure. 
Frequently  the  devices  exhibited  an  open  circuit  at  low  temperatures.  Another  problem 
was  that  data  for  heat  fluxes  higher  than  20  V/cm2  were  unobtainable,  due  to  current 
limitations  in  the  power  supplies. 

Another  concern  was  an  unusual  temperature  response  prior  to  true  nucleate  boiling, 
which  was  not  supported  by  known  conduction  behavior.  Ve  speculated  that  the  response 
might  be  attributed  to  substrate  effects,  poor  convection  at  the  heater  surface,  or  an 
inelastic  interaction  between  the  heater  surface  and  substrate.  Additional  data  now 
indicates  that  the  pattern  of  the  Pt  foil  encourages  the  early  transient  behavior.  A 
detailed  discussion  of  these  results  is  given  later  in  this  report. 

Thick  film  platinum  resistance  devices  on  alumina.  It  was  assumed  that  current 
limitations  would  not  be  encountered  when  using  these  small,  thick  film  devices. 
These  commercially  manufactured  temperature  detectors  have  a  meandering  pattern  of 
Pt,  which  appears  to  be  sputtered  on  the  alumina  support.  The  pattern  could  be 
partially  shunted  to  change  the  overall  resistance  of  the  device  and  thereby  match  the 
desired  impedances.  Ve  obtained  a  resistance  of  10  to  12  ft  at  77  K,  and  the  resistivity 
ratio  Rrt/RRj73  was  0.02.  During  calibration  in  LN2  significant  self-heating  was 
found.  Also,  unusually  small  changes  in  resistance  in  response  to  large  power  inputs 
were  observed.  Calculations  confirmed  the  experimental  results  and  showed  that  907.  of 
the  heat  input  was  absorbed  by  the  alumina.  Several  devices  were  burned  out  during 
preliminary  testing. 

Unsupported  platinum  foil.  Calculations  show  that  a  quartz  or  fiberglass  epoxy 
substrate  could  absorb  up  to  307.  of  the  heat  input  in  liquid  nitrogen  and  107.  in  liquid 
hydrogen.  These  calculations  combined  with  the  data  indicating  the  previously 
described,  unusual  transient  behavior  prompted  the  development  of  an  unsupported 
platinum  foil  heater/thermometer.  Our  experience  with  the  multichannel,  supported 
heater/thermometer  indicated  that  a  single  ribbon  of  platinum  could  be  used  as  the 
test  section.  A  schematic  of  this  configuration  is  shown  in  Figure  5.  Ve  used  various 
devices  of  this  type  to  study  premature  film  boiling  as  a  function  of  size  and  heat 
input.  A  detailed  discussion  of  the  results  is  provided  in  Section  4. 

Carbon  films.  Carbon  has  a  negative  temperature  coefficient  of  resistance,  the 
opposite  of  platinum,  and  a  relatively  high  resistance.  These  attributes  made  it 
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difficult  to  choose  a  carbon  test  section  for  use  with  the  existing  electronics.  For 
example,  the  resistance  of  an  unsupported  Pt  film  at  20  K  often  was  1  ft  or  less,  while 
the  carbon  film  had  resistance  of  nearly  800  ft.  Several  benefits  were  found  from  using 
the  carbon  devices.  First,  a  large  voltage  could  be  obtained  which  could  overcome 
system  noise.  This  meant  a  voltage  response  which  was  up  to  800  times  larger  than  that 
of  the  unsupported  Pt  foil  was  recorded.  A  second  result  was  the  reduction  in  the 
current  required.  The  current  required  was  typically  0.1  A  or  less,  a  value  well  below 
the  5  A  limit  of  the  power  supply.  The  limiting  factor,  therefore,  was  not  the 
magnitude  of  the  heat  input,  but  whether  the  carbon  uevice  could  physically  withstand 
high  power. 

Several  carbon  films  were  tested.  The  first  "film"  was  a  strand  of  fine  carbon 
fibers  stretched  across  a  quartz  substrate.  Its  indefinite  heat  transfer  area  made 
this  device  of  questionable  value.  The  second  film  was  made  of  pressed  graphite  powder 
on  quartz.  This  device  showed  good  sensitivity  when  calibrated  in  the  temperature 
range  within  ±5  K  of  77  K.  However,  the  voltage  response  prematurely  leveled  out 
during  applications  of  pulses  between  5  and  20  V/cm3.  Based  on  previous  transient  LN? 
studies  [1]  a  transition  to  nucleate  boiling  occurs  at  10-3  to  10*1  s;  by  contrast,  our 
results  show  a  steady  state  voltage  occurring  as  early  as  8  x  10'4  s.  Ve  surmise  that 
the  results  are  related  to  the  quality  of  the  carbon  film.  Since  the  film  was  made  by 
rubbing  graphite  onto  quartz,  the  particles  are  only  loosely  held  together  in  a  film  of 
varying  thickness.  Perhaps  this  combination  of  roughness  and  agglomeration  enables 
the  formation  of  "channels"  which  divert  the  flow  of  electrons  through  a  path  of  least 
resistance,  which  decreases  as  more  heat  is  applied. 

The  third  device,  shown  in  Figure  6,  is  a  thin  carbon  film  deposited  on  quartz  by  an 
electron  beam.  The  device  has  been  used  in  previous  heat  transfer  experiments  [7] . 
Its  dimensions  are  approximately  4.11  mm  x  6.35  mm.  Preliminary  results  with  this 
heater/thermometer  submerged  in  LN2  confirmed  earlier  results  reported  by  Steward 
[1] ,  which  warranted  its  use  in  LH2. 

2.4  Data  Acquisition 

A  two-channel  digital  oscilloscope  was  used  to  record  the  voltage  response  of  the 
heater/thermometer  and  the  standard  resistance  wire  from  which  current  was 
determined.  The  oscilloscope  has  ranges  of  200  mV  to  80  V  full  scale.  A  total  of  7964 
points  can  be  obtained  for  each  channel  with  12— bit  resolution  at  100  ns/point.  The 
voltage  responses  were  displayed  on  the  display  screen  and  recorded  on  diskette.  An 
80286-based  desktop  computer  was  used  for  data  analysis  following  each  run.  Data  were 
transmitted  by  RS-232  cable  by  using  a  commercially  available  software  package. 

2.5  Experimental  Procedure 

The  experimental  dewar  was  filled  according  to  the  NIST  guidelines  [6] .  The  dewar 
was  purged  with  gaseous  nitrogen  followed  by  gaseous  hydrogen,  prior  to  filling  with 
liquid  hydrogen.  Safety  procedures  included  posting  warning  signs  donning 
conducting  ankle  straps,  and  arming  the  hydrogen  leak  detection  alarms. 

Electrical  connections  were  made  through  an  air-tight  feedthrough  at  the  top  of  the 
dewar.  A  schematic  circuit  diagram  is  shown  in  Figure  2.  The  test  procedure  began  by 
applying  a  low  voltage  to  the  test  section  until  the  desired  heat  flux  was  verified  by  a 
spot  computation.  The  step  inputs  were  initially  of  10*4  s  duration,  and  incrementally 
increased  up  to  10  s.  The  power  was  likewise  modified  in  small  increments,  and  each 
response  was  recorded. 
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1  igure  6.  Schematic  of  carbon  thin  film  deposited  on  a  quartz  substrate  used  as  a 
h>  ater/thermometer 
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3.  ANALYSIS  OF  TRANSIENT  CONDUCTION  HEAT  TRANSFER 

Immediately  following  a  step  heat  input,  until  the  beginning  of  natural  convection 
and  boiling,  the  heat  transfer  from  solid  to  liquid  obeys  the  law  of  conduction. 
Comparison  of  the  measured  temperatures  with  conduction  calculations  is  a  reliable  way 
of  verifying  the  accuracy  of  transient  temperature  measurements,  and  the  point  of 
digression  marks  the  beginning  of  convection  and  boiling.  The  two  configurations,  (1) 
thin  films  of  Pt  or  carbon  deposited  on  a  solid  substrate,  and  (2)  suspended  Pt  ribbons 
without  a  substrate,  require  two  different  solutions  of  the  transient  conduction 
equation. 


3.1  Thin  Films  on  Solid  Substrates 

Vapor-deposited  films  have  so  little  mass  compared  to  the  substrate  that  the 
sensible  heat  of  the  film,  its  thickness,  and  any  temperature  difference  across  it  may 
be  neglected.  The  following  expression,  for  which  the  variables  are  defined  in  the 
Nomenclature  section  at  the  end  of  the  paper, 


T  = 


ST 


i/z_ 

VMT' 


(1) 


for  the  heater  surface  temperature  of  a  thin  film  was  derived  previously  [1]  and  is 
presented  in  Appendix  A  of  this  report . 

3.2  Ribbon  Heater/Thermometer  without  a  Substrate 

Free-standing  ribbons  are  thicker  than  vapor  deposited  films,  and  the  heat 
capacity  of  the  solid  as  well  as  possible  temperature  gradients  must  be  accounted  for. 
The  two  required  differential  equations  are  for  the  deposited  film  and  the  liquid, 
rather  than  its  substrate  and  the  liquid.  The  following  closed-form  solution  for  the 
temperature  may  be  obtained  through  use  of  Laplace  transforms: 


T' 


+  erfc 


2n+l 


+  z' 


d0  . 


(2) 


A  full  derivation  is  presented  in  Appendix  B. 
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4.  RESULTS  AND  DISCUSSION 


From  measurements  of  voltage  across  the  test  section  of  a  heater/thermoraeter  and 
voltage  across  a  standard  resistor  in  series,  the  power  per  unit  surface  area,  Q,  is 


Q=(Vt)(Vs)/(JtsA), 
and  the  test  section  resistance  is 
h  =  (Vt)(Rs)/Vs  . 

The  test  section  temperature  T  is  found  from 


(3) 

(4) 


T  =  T(Rt),  (5) 
where  T(Rt)  is  the  calibration  function  for  the  particular  heater/thermometer. 

4.1  Transient  Heat  Transfer  to  Liquid  Hydrogen  Using  Carbon  Films 

Calculations  of  conduction  heat  transfer  compared  to  the  experimental  data  carbon 
films  are  shown  in  Figures  7  and  8  for  powers  from  1.15  to  42  V/cm2,  with  heater 
temperature  rise  as  a  function  of  time.  The  data  agree  with  calculations  within 
experimental  error  except  for  a  slight  deviation  toward  the  end  of  the  period, 
possibly  due  to  early  nucleation  at  preferred  sites.  The  deviation  in  slope  of  the 
data  from  the  conduction  calculations  indicates  that  steady  nucleate  boiling  has  been 
reached.  At  7.8  V/cm2,  after  a  period  of  steady  temperature  the  data  exhibit  a  small 
overshoot  at  0.05  s.  Since  the  slope  of  the  temperature  data  returns  to  the  previous 
level,  this  change  in  slope  may  indicate  localized  film  boiling.  The  heat  flux  of  7.8 
V/cm2  corresponds  well  to  the  steady-state  Kutateladze  prediction  [8]  of  9  V/cm2  for 
tne  nucleate  boiling  maximum  or  critical  heat  flux  in  liquid  hydrogen.  For  heat  fluxes 
of  10.8  V/cm2  and  higher,  sharp  transitions  into  film  boiling  can  be  discerned  after  a 
period  of  metastable  nucleate  boiling,  shown  by  a  gradual  increase  in  temperature. 

4.2  Premature  Film  Boiling  Using  Pt  Foil  Heater/Thermometers 

Heater/thermometers  of  supported  and  unsupported  Pt  foil  were  used  to  study 
premature  film  boiling  in  LH2.  Premature  film  boiling,  distinguished  by  a  sudden, 
unexpected  transition  to  film  boiling,  has  been  noted  for  LN2  and  LH2  in  transient 
stepped  power  tests  with  thin  Pt  films.  [4,5]  Ve  never  observed  this  behavior  with  the 
carbon  film. 

An  estimate  for  the  time,  tf,  when  a  transition  to  film  boiling  takes  place  at  a 
solid— liquid  interface  can  provides  a  basis  for  evaluating  premature  film  boiling. 
Schmidt  [2]  developed  a  correlation  for  LHe  for  tf,  which  is  given  by 


tf  =  72 


_k_ 
P  c 


L  ]2 

T 


(6) 


Equation  (6)  is  plotted  on  Figure  9  and  has  a  slope  proportional  to  Q*2;  7  is  computed 
from  carbon  film  experimental  data  and  fluid  properties  for  LH2.  Ve  modified  Equation 
(6)  to  reflect  the  fact  that  tf  approaches  infinity  at  the  critical  heat  flux,  and  to 
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Figure  8.  Results  of  transient  heat  transfer  experiments  from  a  carbon  thin  film  to 
LH?  for  heat  fluxes  ranging  froa  13  to  42  V/cm*. 
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ONSET  OF  FILM  BOILING,  s 


Figure  9.  Experimental  data  from  transient  heat  transfer  experiments  using  a  carbon 
film  in  LH2  which  show  the  time  of  transition  to  film  boiling  for  various  applied  heat 
fluxes.  The  data  are  compared  to  Equations  (6)  and  (7)  which  predict  the  time  for  the 
onset  of  film  boiling. 
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coapute  tf  froa  the  time  of  end  of  the  conduction  period  rather  than  from  the  time  of 
the  beginning  of  the  application  of  power: 


tf  =  tc  +  72 


k 

f  L 

P  c 

»  . 

IHIcr 

(7) 


The  carbon  fila  data  points  show  reasonable  agreeaent  with  Equation  (6).  The  slight 
curvature  in  the  slope  of  the  data  points  at  the  higher  heat  fluxes  is  similar  to  the 
curved  line  of  Equation  (7)  in  Figure  9.  The  agreement  may  be  enhanced  by  optimizing 
7. 

Observations  in  liquid  nitrogen.  Sinha  et  al.  [4]  report  that  following  stepped 
power  inputs,  transitions  to  stable  film  boiling  in  LN2  on  horizontal  0.102  mm  Pt  wires 
were  observed  for  heat  fluxes  as  low  as  407.  of  the  maximum  steady— state  heat  flux  for 
fila  boiling.  Tsukamoto,  et  al.  [5]  photographed  premature  film  boiling  resulting 
from  stepped  power  inputs  to  a  horizontal  0.05  mm  Pt  wire  at  507.  of  the  steady-state 
burnout  heat  flux.  Ve  recorded  experiments  in  LN2  on  high-speed  video  which  used  a 
0.38  am  vide,  unsupported  Pt  ribbon.  Ve  observed  the  development  of  periodic  bubbles, 
which  coalesced  and  enveloped  the  ribbon,  as  well  as  the  corresponding  transitions  to 
film  boiling.  These  transitions  are  unusual  because  they  occur  at  heat  fluxes  from  2 
to  10  V/cm2,  which  normally  support  steady  nucleate  boiling.  [1]  Ve  assume  similar 
behavior  for  liquid  hydrogen,  due  to  the  similarity  in  results. 

Liquid  hydrogen  results.  Ve  found  that  early  and  unpredictable  transitions  to 
film  boiling  in  LH2  occurred  using  unsupported  Pt  ribbons  with  widths  of  0.178  mm  to 
1.524  mm.  The  data,  shown  in  Figure  10,  do  not  follow  the  Q’2  dependence  predicted  by 
Equation  (61  and  exhibit  widely  varying  transition  times  for  the  same  heat  flux.  The 
heater  vidtn  is  on  the  order  of  the  width  of  the  average  nucleate  boiling  bubble.  An 
average  bubble  size,  between  0.246  to  0.52  mm,  and  time  of  departure  from  the  heater 
surface,  between  5  to  36  ms,  has  been  reported  by  Thome  and  Davey,  [9]  for  mixtures  of 
LN2  and  LAr,  and  by  Bevilogua  et  al.  [10] ,  for  bubble  formation  in  boiling  hydrogen.  In 
the  experiments  of  Sinha  et  al.  [4]  and  of  Tsukamoto  et  al.  [5]  the  Pt  wire  diameters, 
respectively  0.102  mm  and  0.05  ma,  were  smaller  than  the  average  bubble  size. 

Premature  fila  boiling  aight  be  a  aanifestation  of  nonunifora  heating  promoted  by  a 
positive  temperature  coefficient  of  resistance  combined  with  the  small  width  of  the 
heater.  In  particular,  if  the  width  of  a  heater  is  nearly  as  small  as  a  nucleate 
boiling  bubble,  a  small  number  of  bubbles  aight  effectively  insulate  a  segment  of 
heater  and  cause  the  temperature  of  that  segment  to  rise  above  the  average.  The  local 
Joule  heating  and  temperature  would  then  increase  further.  The  absence  of  a  nearby 
heat  sink  or  substrate  would  encourage  this  process  to  continue. 

Hvdrodvnaaic  instability.  Ve  apply  the  Taylor-Helmholtz  analysis  to  our  data. 
A  harmonic,  surface-tension  wave  on  the  vapor— liquid  interface  with  a  small  amplitude 
perturbation  may  be  described  by  [11] 

z  =  z0ebt  cos(ax).  (8) 

The  stability  of  the  system  is  determined  by  b.  There  is  a  dominant  wavelength  which 
corresponds  to  the  aaxiaua  value  of  b.  A  maximum  value  of  b  describes  when  the  minimum 
thickness  of  vapor  required  for  stable  film  boiling  has  been  reached  [11] .  The  maximum 
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ONSET  OF  FILM  BOILING,  s 


HEAT  FLUX,  W/cm2 


Figure  10.  Experiaental  data  for  unsupported  Pt  filas  exposed  to  LH2  on  both  sides 
coapared  to  prediction  for  onset  of  f  ila  boiling  presented  in  Equation  6. 
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value  of  b,  44.5  s*1,  corresponds  to  a  dominant  wavelength  of  20  mm.  Comparing  the 
corresponding  time  of  22.5  ms  to  Figure  10  shows  that  many  transitions  occurred  at  even 
earlier  times,  a  reasonable  indication  that  stable  film  boiling  has  been  reached  in 
those  cases.  Other  factors  which  may  account  for  premature  film  boiling  include:  (1) 
reduced  surface  tension  with  a  vertical  test  section,  and  (2)  the  width  of  the  Pt 
ribbon,  many  times  smaller  than  the  dominant  wavelength,  which  encourages  the 
coalescence  of  the  nucleate  boiling  bubbles. 

Temperature  excursions  in  substrate-mounted  Pt  devices.  The  data  of  Figure  11, 
for  substrated-mounted  Pt  devices  in  LHj  differ  from  those  of  the  carbon  film  and  from 
the  unsupported  Pt  foil  in  LH2  in  two  respects:  (l)a  large  temperature  overshoot  at 
the  onset  of  nucleate  boiling,  and  (2)a  brief  temperature  excursion  preceding  the 
overshoot.  These  phenomena  are  associated  with  the  use  of  narrowest  (0.127  mm  wide) 
serpentine  pattern  Pt  heater  mounted  on  a  quartz  substrate.  This  heater  also  had  a 
number  of  imperfections,  such  as  pinholes  and  ragged  cut  edges,  which  could  promote 
the  formation  of  a  hot  spot  as  discussed  earlier.  The  temperature  excursion  can  be 
explained  by  the  development  of  such  a  film  boiling  spot,  and  the  substrate  acted  as  a 
poor  heat  sink  imperfectly  connected  to  the  foil  through  the  epoxy  adhesive.  The 
relatively  large  temperature  overshoot  could  also  be  related  to  this  poor  thermal 
coupling. 


4.3  Estimation  of  Heater  Temperature  Rise  During  Film  Boiling 

For  heat  fluxes  greater  than  the  nucleate  boiling  maximum,  a  period  of 
extraordinarily  low  heat  transfer  resistance,  metastable  nucleate  boiling  ends  with  ?. 
rise  in  temperature  toward  the  film  boiling  level.  For  practical  applications  it  can 
be  assumed  that  the  heat  transfer  coefficient  falls  immediately  to  a  value  of  30  to  50 
mV/(cm2-K),  corresponding  to  film  boiling  [8],  and  the  thermal  properties  of  solid 
materials  control  temperature  rise.  Heaters  with  very  small  mass,  such  as  the  ?t 
ribbons,  jump  quite  abruptly  to  steady-state  film  boiling  temperatures. 

The  solid  lines  for  heat  fluxes  of  13  and  42  V/cm2  in  Figure  12  show  the  results  of 
our  calculation,  using  Equation  (1),  of  temperature  rise  during  the  transition  where 
conduction  into  the  quartz  substrate  controls  the  temperature  rise  and  the  effect  of 
the  fluid  is  neglected.  The  curvature  is  due  to  the  large  variation  of  solid  thermal 
properties.  These  conduction  curves  would  extend  upward  indefinitely  and  therefore 
should  be  discontinued  where  they  intersect  the  steady  state  film  boiling  temperature. 

4.4  Transient  Conduction  Heat  Transfer  from  Pt  Ribbons 

Figures  13  - 16  show  temperature  rise  versus  time  for  heat  transfer  to  saturated  LH2 
from  a  Pt  ribbon  1.778  x  10'2  cm  wide,  3.48  cm  long,  and  4  /xm  thick,  exposed  to  the 
liquid  on  both  sides.  Stepped  power  inputs  were  3,  6,  7,  and  9  V/cm2.  Conduction  heat 
transfer  gives  way  to  boiling  heat  transfer  at  about  1.2  ms  at  3  V/cm2  (Figure  13) ,  and 
the  transition  time  recedes  as  the  power  increases  (Figures  14  -  16).  At  9  V/cm2  the 
transition  to  boiling  begins  at  about  0.4  ms.  In  order  to  describe  the  conduction 
period  realistically,  the  computations  used  temperature  dependent  liquid  and  solid 
properties,  and  the  time  dependent  power  as  determined  experimentally.  The  resulting 
computations,  shown  as  the  solid  lines  in  Figures  13  -  16,  agree  well  with 
measurements  made  during  periods  when  the  power  is  not  constant;  the  variations  in 
power  are  reflected  in  the  different  slopes  at  early  times. 
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Figure  12.  Comparison  of  data  to  calculation,  based  on  Equation  1 ,  of  temperature  rise 
during  transition  to  film  boiling  of  a  substrate-mounted  carbon  film  in  LH2 
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Figure  13.  Comparison  of  experimental  data  for  unsupported  Pt  ribbon  to  computation 
with  temperature  dependent  solid  and  liquid  properties  and  the  time  dependent  power 
for  power  input  of  3  V/cnJ. 


Figure  14.  Comparison  of  experimental  data  for  unsupported  Pt  ribbon  to  computation 
with  temperature  dependent  solid  and  liquid  properties  and  the  tine  dependent  power 
for  power  input  of  6  V/cml. 
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Figure  15.  Comparison  of  experimental  data  for  unsupported  Pt  ribbon  to  computation 
with  temperature  dependent  solid  and  liquid  properties  and  the  time  dependent  power 
for  power  input  of  7  V/cmJ. 


Figure  16.  Comparison  of  experimental  data  for  unsupported  Pt  ribbon  to  computation 
with  temperature  dependent  solid  and  liquid  properties  and  the  time  dependent  power 
for  power  input  of  9  V/cmJ. 
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5.  SUMMARY  AND  CONCLUSIONS 


Ve  studied  transient  heat  transfer  to  carbon  fila  and  Pt  foil  devices  submerged  in 
LH2,  with  emphasis  placed  on  the  onset  of  conduction  and  of  fila  boiling.  The 
experimental  results  are  summarized  in  Table  1  for  each  heater/thermometer  used.  A 
solution  of  Equation  (1)  for  the  transient  conduction  differential  equation  for  a 
massless  heater  bounded  by  a  solid  substrate  on  one  side  and  liquid  on  the  other  was 
used  to  analyze  the  carbon  film  data.  Representation  of  the  unsupported  ribbon  data 
required  a  new  solution,  which  is  presented  in  this  paper  as  Equation  (2).  Except  in 
cases  of  premature  film  boiling,  the  empirical  expression  of  Schmidt  [2],  Equation 
(6),  and  the  modification  in  Equation  (7)  predict  tf  at  the  higher  heat  fluxes. 
Premature  film  boiling  can  be  related  to  the  narrowness  1.524  mm  or  less  of  the  heaters 
relative  to  the  size  of  nucleate  boiling  bubbles,  irregularities  in  the  heater 
material,  and  a  positive  coefficient  of  electrical  resistance. 


Table  1 .  Summary  of  Transient  Heat  Transfer  to  Liquid  Hydrogen 


Type 

of 

Heater 

Transient 

Conduction 

Period 

Equation 
for  Trans. 
Conduction 

Time  to 

Onset  of 

Film 

Boiling 

Equation 
for  Onset 
of  Film 
Boiling 

Comments 

or 

Exceptions 

Carbon 

Thin 

Film 

10  /xs  to 

20  ms 

Eq(l) 

10  ms  for 
10.8  V/cm2 

Eqs  (6) 
and  (7) 

Critical 
Heat  Flux 
7.8  V/cm2 

Pt  Foil 
Ribbons 
Boiling 

10  /is  to 

1  ms 

Eq  (2) 

3  ms 
for 

7  V/cm2 

Eqs  (6) 
and  (7) 

Premature 

Film 

related  to 

narrow 

heaters 

Pt  Foil 
on 

Quartz 

1  (iS  to 

5  ms 

Eq(l) 

20  ms 
for 

12  V/cm2 

(not 

studied) 

Temp. 

excursions 
due  to 
hot  spots, 
substrate 
effects 
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7.  NOMENCLATURE 


A  surface  area,  cm 

b  growth  parameter 

c  specific  heat,  J/(g  K) 

k  thermal  conductivity,  W/(cm  K) 
L  heat  of  vaporization,  J/g 

m  wave  number  =  2x/A 

Q  applied  heat  flux,  W/(cm3  s) 

Qcr  critical  heat  flux,  W/(cmJ  s } 
q  heat  flux  per  unit  volume,  W/cm3 
Rs  resistance  of  standard  resistor,  fi 
Rt  resistance  of  test  resistor,  ft 

t  time  s 

t'  dimensionless  time  =  taJ6i 

tc  time  of  conduction  period,  8 

tf  time  of  onset  of  film  boiling,  s 
T  heater  surface  temperature,  K 

v,i  denote  vapor,  liquid 

8,f  denote  solid,  fluid 

erfc 


T'  dimensionless  temperature 

=  Tas/(q£3) 

Vt  voltage  across  the  test  section,  \ 

Vs  voltage  across  the  standard 
resistor,  V 

x  length  of  fluid— vapor  interface,  cm 

z  amplitude  of  perturbation,  cm 

zo  perturbation  constant 

a  thermal  diffusivity  =  k/(p.cl,  cin*/s 
0  dimensionless  parameter 

=  (ks/kf)(as/af) 

6  ribbon  thickness,  cm 

7  fit  factor 

p  density,  g/cm3 

A  wavelength,  cm 

a  surface  tension,  dyn/cm 

0  integration  variable  for  time 


complementary  error  function 
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Appendix  A 


TRANSIENT  CONDUCTON  HEAT  TRANSFER  - 
FLAT  PLATE  HEATER  ON  A  SOLID  SUBSTRATE  WITH  LIQUID  COOLANT 


Assume: 

1.  Infinite  flat  plate  heater  in  the  x  —  y  plane  at  z  =  0.  Note:  If  the  substrate 
were  finite  in  the  — z  direction  one  would  have  to  make  an  assumption  about  heat 
transfer  at  a  bottom  boundary.  The  results  are  the  same  until  a  temperature  wave 
penetrates  the  finite  substrate. 

2 .  Heater  has  zero  thickness . 

3.  Stagnant  liquid  hydrogen  lies  above  the  plate  (z  >  0);  no  forced  convection  is 
present . 

4.  A  solid  substrate  lies  below  the  plate  (z  <  0) . 

5.  Times  of  concern  are  short  enough  that  thermal  convection  currents  are  not 
established  and  no  bubbles  have  formed.  This  would  mean  non— equilibrium  of  phases 
(superheating)  in  cases  where  Tf  >  T8aturation* 

6.  Heating  is  at  constant  volume,  i.e.  no  work  of  compression. 

7.  Radiation  heat  transfer  is  much  less  than  conduction  heat  transfer. 

8 .  Fluid  and  sol  id  propert  ies  are  constant . 


The  Fourier  heat  conduction  equation  in  1  dimension  is  presented  for  the 
following  cases: 

A.  Solid  Substrate  (subscript  8) 


dig  1  dig 

tt^a-gW  z-0’  (1) 


where 

T* 

=  temperature  of  substrate  at  z  and  time  6 
=temperature  before  disturbance. 

0 

=  time 

Ks 

Os 

=  thermal  diffusivity  =  — r — 

P  stps 

=  thermal  conductivity 

Pa 

=  density 

Cp8 

=  specific  heat 

Boundary  Conditions  for  the  Substrate 

T,(M)  =0 


(2) 
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(3) 


lim  Ts(z,0)=O 


P&  -  Kg 


the  part  of  the  heater 
pover  being  dissipated 
in  the  substrate. 


B.  Conduction  Equation  for  the  Fluid 


(4) 


d2Tf  1  5Tf 
dz*  ~  Of  Iff) 


where  Tf  =  temperature  of  fluid  at  z, 

0  =  tenperature  before  disturbance. 

Boundary  Conditions  for  the  Fluid 


Tf(z,0)  =0 


lim  Tf(z,0)  =  0 

Z-KB 

ran 


Pf=  -if 


0,0 


part  of  the  heater  power 
being  dissipated  in  the 
fluid. 


C.  Condition  at  the  substrate/heater/liquid  interface: 


total  heater  power  — 
P  =  Pf  +  Ps  known  parameter 


T,(<M)  =Tf(O,0) 
D.  Solution  for  the  Substrate 


(5) 


(6) 

(7) 

(8) 


(9) 

(10) 


Laplace  transform  with  respect  to  0 
d2t  1 

W  (z’s)  =  q7  [st(z’s)  — T(z,°)]  .  (11) 

FromB.C.  #2T(z,0)  -*0, 

where  t(z,s)  =  L[  T(z,0)] 

s  =  the  Laplace  transform  parameter. 

The  solution  of  (11)  is 


t  =  Ci  ez'/®7<*s  +  Cj  e-*/i75  . 
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FromB.C.  #3  Ca  =  0.  Therefore, 


t  = 

Transfora  B.C.  #4: 


(12) 


P* 

r 


tdt 

ft 


1 0 ,  s 

Then  from  Equation  (12) , 

P8 


~  =  KgCi  <Jafch  Qz^*ft**  z=0 


or 


P»1  r-r- 
Cl  “  5"  kT  . 


Therefore, 


t=i 


P« 


yfclse 


zjsjcta 


(13) 


It  is  assumed  in  the  following  development  of  equations  that  Ps  is  constant. 
Further  explanation  is  provided  in  the  addendum  following  this  appendix. 


Use  the  following  transform  relations: 


L 
and 


tz3/4Qg0 


a 


(14) 


(15) 


Then  combining  Equations  (13)  and  (14)  yields  the  following  expression: 


^s1 


p. 
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From  Equation  (15) 


t  =  L 


-6  p 

f-  Va»/(*T)  dr 

0ls 


(where  t  =  dumay  variable) . 


Therefore, 


J 

e-z3/(4ogr) 

0  H 


dr  . 


(16) 


A  change  of  variables  is  now  introduced: 


Let  A  = 


2V5s r 

Therefore  -i-  =  }  and  r  =  — — 

v'r  z  4A3Og 

Hence,  dr=  =  -±*1Q 

4  A  a*  2  A3  o8 


Also, 


r  =  0  =»  A  = 
r=  0*  A  = 


®  » 
z 


2jc£0 

The  following  expression  can  be  derived: 

z 

_ I  2Vag0 


Tg  =  V^s/* 


— AJ 


2y/5g 


1  Z3 

r?r*j 


dA  , 


or 


Tg  =  ^4-^ 
/* 


dA  . 


2^0*0 


(17) 
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Now  we  perform  an  integration  by  parts:  Judv  =  uv  -  Jvdu. 


Let  dv  —  ^  and  v  =  — ^ 


u  =  e  ^  and  du  =  -2Ae  ^  dA 


Then, 


Ta  =  h^. 
I.V5F 


f£] 


r® 


or 


P» 
T»  = 

*8 


-Z  ' 


0D 

-2 

z 

e~; 

z 

2  Va*0 

2jc£0 

■0D 

-9  2z 

~  7* 

e- 

z 

2jaj) 

dA 


z-erfc 


IsfaJ) 


Therefore, 


Ts  =  i4 


K, 


2yfc6Ti<z*%IAa*e- z-erfc 


2y]az(r 


at  z  <  0  . 

Substrate  Surface  Temperature  (z  =  0) 

T,(O,0)  =2  h&Jd. 
ft  *8 

E.  Solution  for  the  Fluid 


(18) 


(19) 


By  starting  with  eq.  (5)  and  boundary  conditions  (6),  (7),  and  (8)  we  get  the 
same  solution  as  (18)  except  the  properties  are  fluid  properties  and  the  power  is  Pf . 


Pf 

Tf  =  f? 


2Vo£*7?e~*3/4af0_  z-erfc 


2/Etf9j 
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(20) 


and  at  the  liquid  surface  the  temperature  is: 


2Pf  <Ja 7 

Tf(°>^)  =  75TTf  ^  •  (21) 


F.  Solution  in  terns  of  P  and  0 

At  the  interface  Ps  +  Pf  =  P,  and  Tg(O,0)  =  Tf(O,0)  =  T(O,0).  Therefore,  from 
Equations  (19) ,  and  (21) : 


P  =  I(o1ei^ 


Jf' 

V5f. 


(21a) 


ADDENDUM  -  PROOF  THAT  Pg  AND  Pf  ARE  CONSTANT 

From  Equation  (12)  we  know  that  t  =  Cj  e2^8/**8  for  the  substrate  (z  <  0) . 

For  the  fluid  t  =  Cje-2^8/0^  ,  (z  >  0)  .  (A) 

At  z  =  0  the  temperature  of  the  substrate  =  tenperature  of  liquid,  therfore  Ci  =  C3  =  C. 

P8  =  s  C  IfVs/Of  , 
and 

Pf  =  s  C  Igv/s/ag  . 

Therefore, 

Pg  If  ag 

FT  =  =  (co"8t)i. 

P  =  Pf  +  Pg  =  (const) j  (given)  ,  (B) 

P  =  Pf ( 1  +  const i)  =  (con8t)2  . 
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(const) j 

Siailarily  Pg  =  (const<) . 


It  is  not  necessary  to  know  that  P8  is  constant.  Using  condition  (A)  above  and 
Equation  (22)  gives  the  saae  result. 


Appendix  B 


TRANSIENT  CONDUCTION  HEAT  TRANSFER  FROM  A  RIBBON  HEATER 
VITH  INTERNAL  HEAT  GENERATION 


The  unsteady  state  conduction  equation  is 

q  1  dl 
«*• 

where  q  =  generated  heat  flux/volume 

k  =  thermal  conductivity 
o  =  thermal  diffusivity 
T  =  temperature  -  initial  temperature 
t  =  tine 


1) 


A.  Development  of  equations  for  the  solid 
Equation  (1)  for  the  solid  is: 

^T  q  1  5T 

3yi  +  I7  =  573t-  (2) 

Boundary  Conditions: 

T(y,0)  =  0  (initial  condition).  (3) 

dl 

3t(0,t)=0  (symmetry). 

Change  to  dimensional  ess  variables: 

y 

Let  y'  =  ^  (dimensionless  space  variable) . 

To  write  Equation  (2)  in  dimensionless  terms,  the  following  must  exist: 
dn  q  T 

W  =  ^  =  y7' 

T 

Let  T'  =  j  . 

T'A  q  q<f» 

Then  (y'S)3  t=J  fj  or  A  =  IT"  • 


31 


1 


t 

Let  t'  =  g  .  Then, 

1  dl  1  AT'  q 

H  aw  H  rr- 

T'  q^ 

Since  p-  [=]  1  and  7  =  | —  , 

1  qtf*  1  q 
al7*B  =  K7  ’ 

where  B  =  — ^  . 

a 

Recap:  The  dimensionless  space  coordinate,  temperature,  and  time  are: 


y'  =  3  or  y  =  fiy1 
TIS 


T'=*75  or  T  =  T,fr 


t'  = 


tag 

W 


q  6* 

IT 

t'6* 


or  t  = 


a« 


(4) 

(4) 

(4) 


When  these  variables  are  substituted  into  Equation  (2)  we  have  the  dimensionless 
equation  for  the  solid: 


dJT'  ffl' 

3r*+1  =  W 


(5) 


B.  Development  of  equations  for  the  liquid 

The  conduction  equation  for  the  liquid  lacks  the  heat  generation  term 
d*T  1  dl 

W*  ’  01 5t  .  <6> 

The  dimens  ionless  equation  for  the  liquid  is 
d*V  a ,  dl' 

W*  =  3TJE7-  <7> 

Boundary  condition  for  the  liquid 

T(®,t)  =  0  .  (8) 
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C.  Coupling  the  equations 

To  couple  Equations  (2)  and  (6) : 

at  y  =  6/2  TB  6/2,t  =Ti$/2,t  , 

(9) 

fdT  r  m 

and  Ig  7j-  =  Ki  75— 

^y  =6/2  l>^Jy  =  «/2 

(10) 

Or,  in  terms  of  the  dimensionless  variables: 

y' =1/2  V(l/2,t')=Tl'(l/2,t') 

(9') 

and'-[f:]y,=1/2  =  I‘(f:]y,=i/2 

(10') 

D.  Transformed  Solution 
Laplace  Transform 


53r  1 

Solid:  3pT  +  i  =  sr-  r(y/  ,0) 

Initial  condition.  r(y/  ,0)  =  0) 
d>r'  1 

W*~8T  =  ~* 

where  r  is  the  Laplace-transformed  temperature 

s  is  the  Laplace  transform  variable 

1 

Particular  Solution  rP  =  p  • 

Homogeneous  Solution  rc  =  Cje^y  +  C4e~v^y  . 

Total  solution  r  =  Cse^®y  +  C4e"~^®y  +  s3  •  (11) 

dr  ' 

The  symmetry  B.C.  -*  -jr-r 

L  °y  Jy'  =o 

From  Equation  (11) ,  ^  ^  =  [  Cj^se^'  -  C4Vse^®y'  ]y/=Q  , 

which  gives 
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Cj  =  C4  . 


Therefore, 


d2T  a g  f  yO] 

For  the  liquid  -^n  =  ”  (sr-  r(y'<0)J 

r  =  Cott*  ^  , 


r(a,t)  =  0  or  t(od,s)  =0  — *  Ci  =  0  . 

Therefore , 

T.C*-«r'V575.  (13) 

C3  and  Cj  can  be  evaluated  from  Equations  (12)  and  (13)  and  the  coupling 
condition  from  Equation  (9')  at  y'  =  1/2: 

[e*£/2  +  e~^®/2j  +  p  =  fte"^2^5*  .  (14) 

To  use  Equation  (10') ,  differentiate  Equations  (12)  and  (13)  and  evaluate  at 
y' =1/2: 

C3  (e^/2  +  e"^/2]  =  -C^  e-Vs/2^7^  5 

C,  ^  VS75T  [e^/2  +  e-^®/2]  =  CJe~^/2'/^  .  (15) 

By  combining  Equations  (14)  and  (15)  we  obtain  the  following  expression: 

cje’fr'2  ♦  e^5/2]  ♦  ^  =  V57 of  C  ,[*W2  *  e^/2] 


eVs/2  +  e-Vs/2j  +  ^fi/2  _  e^/2 


It  Qi  (KpCp) g 

Aere^^j—  = 


34 


By  combining  Equations(12)  and  (16)  the  following  expression  is  obtained: 

1  r  e^'  ,  »-&'  j 


a’  [  US/*  +  e-Vs/2]  _  e^5/2 


2Vsy'  +  g^/sy' 


s7  |  Ae^/2  _  Be^®/2 

where  A  =  0  +  1 

B  =  /?—  1. 

By  long  division 


Lf* 

p  Lx 
■ 


1  e/s(y'-l/2)  +  B  eVs(y'— 3/2)  +  B’  eVff(y#-5/2)  + 

x~  p  p 


i  P-^(y,+1/2)  +  B  e— Vs(y/+3/2)  +  B2  e~Vs(y/+5/2)  + 

X  P  P 


1  _nfe-^(1/2-y/)  +  Be^3/2-y,5 
P  X  S  L  i  I  8 

fBl2 e~^^^2_y  ^  e~/s(l/2+y  ) 


+  .  .  .  + 


B  e-v®(3/2*y')  rB1Je-««(5/2-y') 

+  +1  +  • 
I -  X  s 


1.  Tiae  Doaain  Solution 

The  inverse  transform  for  the  solid  is: 


T'  =  t'  - 


f*  [  erfc 

1/2  - 

*  B  erfc  3/2  -  y'} 

J0  [ 

[“279 

i  ~2JV 

fBl 2  , 

+  erfc 

I 


+  ®  erfc 


l-r 

2 


+  .  .  .  +  erfc 


1/2  +  y' 

“275 


'3/2  +  y/- 

+ 

B 

2 

erfc 

s/2  ♦  yr 

+  *  .  • 

[“275 

I 

Ub  \ 

d  e 


or 


T'  =  t'  - 

+  erfc 


'~ii  m-rfeHcf— 

I  n=0  |lj  JO  [  [  275 


(17) 
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